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PREFACE 


The National Ecology Center of the U.S. Fish and Wiidlife Service (FWS) 
is supporting a series of field research studies to document relationships 
between hydric soils and wetland vegetation in selected wetlands throughout 
the United States. This study is one of that series. It is a continuation of 
the FWS effort, begun by Wentworth and Johnson (1986), to develop a procedure 
using vegetation to designate wetlands based on the indicator status of wetland 
vegetation as described by the FWS "National List of Plants that Occur in 
Wetlands" (Reeu 1986b). This list classifies vascular plants into one of five 
categories according to their frequency of occurrence in wetlands. Concurrent 
with the development of the wetland plant list, the Soil Conservation Service 
(SCS) developed a National list of hydric soils (SCS 1985b). Studies supported 
by the National Ecology Center quantitatively compare associations of plant 
species, designated according to their hydric nature using the Wentworth and 
Johnson (1986) procedure, with the hydric nature of soils according to their 
designation on the SCS hydric soils list. The studies are being conducted 
across moisture gradients at a variety of wetland sites throughout the U.S. 
Several studies have been modified to obtain concommitant information on 
groundwater hydrology. 


These studies were conceived in 1984 and implemented in 1985 in response 
to internal planning efforts of the FWS. They parallel, to some extent, 
ongoing efforts by the SCS to delineate wetlands for Section 1221 of the Food 
Security Act of 1985 (the swampbuster provision). The SCS and FWS provided 
joint guidance and direction in the development of the Wentworth and Johnson 
(1986) procedure, and the SCS is currently testing a procedure that combines 
hydric soils and the Wentworth and Johnson procedure for practical wetland 
delineation. The efforts of both agencies are complimentary and are being 
conducted in close cooperation. 


The primary objectives of these studies are _ to: (1) assemble a 
quantitative data base of wetland plant community dominance and codominance 
for determining the relationship between wetland plants and hydric soils; 
(2) test various delireation algorithms based on the indicator status of 
plants against independent measures of hydric character, primarily hydric 
soils; and (3) test, in some instances, the correlation with groundwater 
hydrology. The results of these studies also can be used, with little or no 
supplementary hydrologic information, to compare wetland delineation methods 
of the Corps of Engineers (1987) and the Environmental Protection Agency 
(Sipple 1987). 


Any questions or suggestions regarding these studies should be directed 
to: Charles Segelquist, 2627 Redwing Road, Creekside One Building, Fort 
Collins, CO 80526-2899; FTS 323-5384 or Commercial (303) 226-9384. 
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INTRODUCTION 


Natural landscapes are becoming increasingly difficult to preserve in the 
United States. Even ecosystems exposed to some disturbance and/or manipulation 
are limited and under pressure for urban, industrial, or agricultural develop- 
ment. Wetland systems are among those receiving the greatest pressures. 
Crumpacker (1984) estimates that 70%-90% of the natural riparian ecosystems in 
the United States have been lost to human activities. As a result it has 
become necessary for Federal, State, and private organizations concerned with 
acquisition and management of wetland systems to devise efficient and accurate 
techniques for the assessment of biotic and abiotic characteristics of wetland 
Sites. Activities along these lines have greatly increased over the past 
several years. One example is a volume edited by Warner and Hendrix (1984) on 
California riparian systems ecology, which includes relationships of wetland 
plants with soils and abiotic features. 


The U.S. Fish and Wildlife Service is one Federal agency actively pursuing 
the refinement of inventory technologies and methodologies for classifications 
and data analyses. The Service defines wetlands as follows (Cowardin et al. 
1979): 


. transitional between terrestrial and aquatic systems where the 
water table is usually at or near the surface or the land is covered 
by shallow water . . . wetlands must have one or more of the follow- 
ing three attributes: (1) at least periodically, the land supports 
predominantly hydrophytes; (2) the substrate is predominantly 
undrained hydric soil; and (3) the substrate is nonsoil and is 
Saturated with water or covered by shallow water at some time during 
the growing season of each year ... The upland limit of a wetland 
is designated as: (1) the boundary between land with predominantly 
mesophytic and xerophytic cover; (2) the boundary between soil that 
is predominantly hydric and soil that is predominantly nonhydric; or 
(3) in the case of wetlands without vegetation or soil, the boundary 
between land that is flooded or saturated at some time each year and 
land that is not. 


Hydric soils are defined by the Soil Conservation Service (1985a) as 
soils that in an undrained condition are saturated, flooded, or inundated long 
enough during the growing season to develop anaerobic conditions that favor 
the growth and regeneration of hydrophytic vegetation. 


Reed (1986a) compiled a list of wetland plants of New Mexico, assigning a 
wetland indicator number to each species based on its frequency of occurrence 
along the moisture gradient. Wentworth and Johnson (1986) devised a system of 
using these wetland indicators for the designation of wetlands. 








In 1985, the U.S. Fish and Wildlife Service initiated studies to determine 
the degree of correlation of wetland plant indicator ratings (Reed 1986b) and 
the SCS hydric soils list, test the validity of the Wentworth-Johnson system, 
and assess the validity of Reed's species ratings. The studies include various 
types of wetland systems found in the United States. Our study was conducted 
on western riverine systems in New Mexico. Soils were classified and 
delineated by Mr. Jimmy M. Gass, Soils Scientist for the U.S. Forest Service 
in Albuquerque, New Mexico. 


DESCRIPTION OF THE STUDY AREAS 


We conducted this study on reaches of the Gila and San Francisco Rivers 
in New Mexico (Figure 1). There have been some studies of riparian vegetation 
on intermittent and perennial streams in this portion of New Mexico and 
adjacent Arizona. Glinsky (1977) described the effects of livestock grazing 
and stream bed erosion on the distribution and regeneration of cottonwoods and 
sycamores along Sonoita Creek in Arizona. Turner (1974) compiled four maps of 
riparian vegetation along the upper Gila River in Arizona and documented 
changes in channel width and vegetation composition, possibly in response to 
dam construction. Hubbard (1977) conducted a biological inventory of the 
lower Gila River valley in New Mexico. Along the Mimbres River in southwestern 
New Mexico, Boles (1978) identified a pattern of community and species replace- 
ment due in part to elevational change. Egbert (1981) conducted a survey of 
the flora and fauna of the Gila Riparian Preserve of The Nature Conservancy. 
Medina (1984) analyzed riparian plant communities and soils of three 
intermittent creeks in the Fort Bayard watershed in New Mexico. The relation- 
ships of soil physical characteristics and pedogenesis to vegetation were 
addressed by Brock (1985) at four sites on the Gila River and two sites on the 
San Francisco River in New Mexico. Community structure of riparian vegetation 
at three sites along the Gila River in New Mexico was analyzed by Hardesty 
(1986). 


Many of the dominant riparian species presently occurring along the Gila 
and San Francisco Rivers are apparently relicts of the deciduous forests 
(Arcto-Tertiary Geoflora) that covered the entire area approximately 15 million 
years before present (Hardesty 1986). Brown (1982) stated that many present- 
day Southwestern riparian trees and shrubs are the same species that have been 
present throughout the Southwest for several million years. Dependence on the 
mesic riparian environment is the common bond uniting the relictual species. 
Henry (1981) hypothesized that the lowering of vegetational zones during the 
Pleistocene had little effect on riparian vegetation in the canyons and smaller 
mountain ranges of the Southwest. However, Hardesty (1986) felt that along 
the Gila River, a definite change in riparian species composition occurs over 
an elevational change of 100 m, and it seems likely that some of the riparian 
vegetation patterns along the Gila may have been affected by the 900-1200 m 
downward displacement of vegetational zones during the Pleistocene, as 
suggested by Antevs (1954). 
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Maker et al. (1978) indicated that most of the soils along these rivers 
in New Mexico are of the Haplargrids-Rough Broken Land type. The topography 
associated with these soils is sloping alluvial fans and terrace tops and can 
be very steep on upland ridges and terraces. 


The general elevation of the study areas is 1500 m + 100m. Water 
temperature in the two streams averajes about 21 °C, and mean annual 
precipitation is 250 mm. There are approximately 190 frost-free days per 
year, and the mean winter air temperature is 5 °C, with a summer mean of 
22 °C. This description has been adapted from Monthly Climatic Summary, 
Office of State Climatologist, New Mexico Department of Agriculture, Las 
Cruces, New Mexico. 


The vegetation of the reaches studied is typically subjected to periodic 
disturbances. These disturbances are mainly of two types: natural severe 
flooding from extensive steep watersheds has occurred at 15- to 20-year 
intervals; and, in places, there still is intensive livestock grazing along 
the rivers. These disturbances have resulted in much of the riverside 
vegetation being in varicus stages of succession. The type and degree of 
disturbance is included in the results section of this report. 


METHODS 


DATA GENERATION 


Four dominant soils were identified along the two rivers by Gass (pers. 
comm.). Soils were not classified to the series level; however, hydric soils 
were identified based on available information at the request of the FWS. The 
study areas were extensive enough to permit four replications and five plots 
per replication on each of the soils, on each of the rivers (Table 1). Pits 
were dug in each soil type and a soil analysis performed. We placed five 
replications rather than four in lower terrace soils on both rivers, in the 
event that some plots in a replication should fall outside the soil class. 
Four ex.ra replications were placed on the swale soils on the Gila River 
because there was extensive marsh area on one portion, possibly indicating a 
different hydric soil. 


Field work was conducted during June and August of 1986. Vegetation was 
sampled on each of the five randomly selected plots within each replication. 
The field sampling technique (Table 2) followed U.S. Fish and Wildlife Service 
guidelines for this study. All plants were identified to species. All plants 
found and their frequency of occurrence index numbers, obtained from Reed 
(1986a), are included in the Appendix. A number of species commonly found on 
the plots were not on the New Mexico list (Reed 1986a). In order to analyze 
and compare all data from all plots, species that were not on the New Mexico 
list were assigned an index number based on their frequency of occurrence in 
this study, associations with species that were on the list, and known riparian 
characteristics from the literature. 





Table 1. Location of replications on soil classes of the Gila and San 
Francisco Rivers. Numbers in parentheses represent the order that 
replications were sampled in the field. 




















Gila River ____ San Francisco River 
Soil Rep. Range Twnshp. Sec. Rep. Range Twnshp. Sec. 
Upper 1(11) R1I7W T17S = 17NE1/4 1(20) R20W 711S 27SE1/4 
terrace 2(12) R17W T17S = 17NE1/4 2(21} R20W T9S 27NE1/4 
(Fluventic 3{13) RI7W T17S 8SE1/4 3(22) R20W T9S 27NE1/4 
Ustochrept) 4(17) R17W~ T17S 8SE1/4 4(23) R20W T9S 27NE1/4 
Lower 1(2) R17W T17S = 28NE1/4 1(18) R20S 7118S 27SE1/4 
terrace 2(5) R17W T17S_ 17NE1/4 1(19) R20W T11S 27SE1/4 
(Typic 3(10) R17W T17W~ 16NW1/4 3(24) R20W T9S 27NE1/4 


stifluvent) 4(14) R16W 7115S 6NW1/4 4(25) R20W 192 27NE1/4 
5(15) Rl6W T15S 6NW1/4 5(34) R20W 1711S 4SE1/4 


Sandbar 1(1) R1I7W  T17S = 28NE1/4 1(35) R20W 1711S 4SE1/4 
(Aquic 2(3) R1i7W T17S ~~ 16NW1/4 2(36) R20W 1711S 4SE1/4 
Ustifluvent) 3(4) R17W  T17S 9Sw1/4 3(37) R20W 1711S 4SE1/4 
4(16) Rl6W 1714S 31SE1/4 4(38) R20W 1711S 4SE1/4 


Swale 1(6) R17W~ T17S 9SW1/4 1(26) R20W T9S 27SW1/4 
(Typic 2(7) R17W~ T17S 9SE1/4 2(27) R20W T9S 27NW1/4 
Fiuvaquent) 3(8) R17W = T17S = 21SW1/4 3(28) R20W T9S 27NE1/4 
4(9) R17W T11S = 21NE1/4 4(29) R20W T9S 27NE1/4 


1(30) R20W 7110S 8SW1/4 
2(31) R20W 110S 8SW1/4 
3(32) R20W T10S 8SW1/4 
4(33) R20W  T10S 8SW1/4 





DATA ANALYSIS 


Weighted averages for each life-form stratum within each plot were 
computed using the Cornell ordination program (Gaugh 1977). The ORDIFLEX 
equation is: 


p 
I; ES /O I. .) 


i=] 'Y 





Table 2. Vegetation sampling design for the Gila and San Franciso Rivers. 





Sample units per 
Stratum name Stratum identifier Measurements rep. N/size 








Trees dbh >7.5cm N individ./species 5/100m? 
dbh individuals 





Tall shrubs dbh <7.5cm N individ./species 5/4m? 
ht. >1.3m N leaders/individ. 
Short shrubs ht. >0.5m to <1.3m N individ./species 5/4m? 
Ground cover ht. <0.5m (woody) % classes 10/0.5m? 
all herbaceous sp. (Daubenmire 1968) 
where W, = weighted average for a plot within replication 
I. = "importance" value for species "i" in plot "j," where the 
J importance value is dbh quantity for tree stratum, density for 
shrub strata, and % cover class for ground cover stratum 
ar = frequency of occurrence index number by species 
p = number of species in a stratum within a plot 


Each stratum was analyzed independently, producing a maximum of four 
separate weighted averages for each plot within a replication. Data were 
analyzed by strata (life-forms), instead of combining strata, because we 
believe that it added important ecological information on differences among 
life-forms, and also better accommodated the different sampling methods used 
for the different strata. 


An analysis of variance (ANOVA) was performed to determine if there was 
an overall significant difference between soil classes within each stratum and 
when all strata were combined. Mean weighted averages of the five plots per 
replication were used here instead of the pooled mean of all plots (independent 
of replication). For intersoil relationships, a Duncan's Multiple Range Test 
(DMRT) was performed to determine specific significant differences among the 
means within a stratum and overall. Both ANOVA and DMRT were computed using 
the General Linear Model (GLM) procedure provided by the Statistical Analysis 
System (SAS) package (Ray 1985). 
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The same procedure was followed to analyze the data using presence/absence 
data instead of importance values. The I. values in the above equation were 


either replaced with a "1" or "0," denoting presence or absence. 


RESULTS 


The floodplain land-forms of the two rivers were described by Gass (pers. 
comm.). A typical profile of these land-forms and the subsurface features is 
depicted in Figure 2. The soils associated with the floodplains of the study 
areas are derived from recent alluvium; most commonly Gila conglomerate. The 
mean annual soil temperature is approximately 13 °C. 


The four floodplain soils, ranging from higher (driest) to lower 
(wettest), are Fluventic Ustochrept, Typic Ustifluvent, Aquic Ustifluvent, and 
Typic Fluvaquent; the last three were classified as hydric soils (Table 3). 
We have assigned nontechnical synonyms to the soils subgroup classes: upper 
terrace soil to Fluventic Ustochrept, lower terrace soil to Typic Ustifluvent, 
sandbar soil to Aquic Ustifiuvent, and swale soil to Typic Fluvaquent. Water 
tables vary among these very deep (200 cm) soils and are reflected in their 
classification. The Typic Fluvaquent is associated with ponded water and an 
indicator of this property is the odor of hydrogen sulfide. The other two 
hydric soils are associated with flowing water. While the three lower soils 
were designated hydric, quantitative seasonal] data on depth to water, soil 
moisture, flooding frequency, and flooding duration were not available. Thus, 
these hvdric designations should be considered tentative and may need to be 
adjusted as more data become availal e. 


The Fluventic Ustochrept is influenced by overland flow of water but is 
stable enough (rarely flocded) to develop a cambic (bw) horizon. Recharge of 
soil water is from direct precipitation plus overland flow. Recharge of soil 
water on the remaining soils is from direct precipitation, overland flow, and 
high water table. 

Following are descriptions of the four soils: 

Fluventic Ustochrept, loamy-skeletal, mixed mesic. 


Pit location: 1T178,R17W,SE1/4,Sec.2.Gila River. 








Horizor 
Designation Description 
0 3 to 0 cm; litter layer of undecomposed leaves, twigs, etc. 
A 0 to 9 cm; dark grayish brown (10YR 4/2) extremely 


gravelly sand loam, very dark grayish brown (10YR 3/2) 
moist; moderate fine granular structure; soft; many fine 
interstitial pores; common very fine roots; 75 percent 
rock fragments; neutral; abrupt smooth boundary. 




















A - Water, current channel 


B - Linear sandbar 


- Erosion slope 


- Gila conglomerate 
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Figure 2. 


Water and saturated 
sediments 


Generalized profile of floodplain on the Gila and San Francisco 


Rivers, including landforms. 
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Bw 


Cl 


9 to 19 cm; grayish brown (10YR 5/2) very stony sandy loam, 
brown (10YR 4/3) moist; weak subangular blocky parting to 
weak fine granular structure; soft; many very fine and fine 
roots; many fine tubular pores; 60 percent rock fragments; 
neutral; abrupt boundary. 


19 to 44 cm; dark grayish brown (10YR 4/2) very stony sandy 
loam, very dark grayish brown (10YR 3/3) moist; moderate 
medium and coarse subangular blocky structure; many very 
fine and fine, common coarse roots; many fine tubular 
pores; 50 percent rock fragments; neutral; abrupt wavy 
boundary. 


44 to 105 cm; brown (10YR 5/3) extremely gravelly sandy 
loam, yellowish brown (10YR 3/4) moist; massive; soft; 
many very fine and fine, common medium and coarse roots; 
many very fine and fine tubular pores; 70 percent rock 
fragments; neutral; abrupt wavy boundary. 
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Table 3. 


Francisco Rivers. 





High water table 











Soil and site characteristics for sample areas on the Gila and San 











Moisture Temp. Depth —__—s«éFlooding 
Soil regime regime (cm) Months Freq. Duration Months 

Fluventic 

Ustochrepts Ustic Mesic 200 coe Rare Very June- 

(upper brief Sept 

terrace) 

Typic 

Ustifluvents Ustic Mesic 150- Jan- Common Brief May- 

( lower 200 Dec June, 

terrace) July- 
Sept 

Aquic 

Ustifluvents Ustic Mesic 50- Jan- Common Long May- 

(sandbar) 150 Dec June, 
July- 
Sept 

Typic Aquic Mesic 50 Jan- Frequent Very May- 

Fluvaquents long Sept 

(swale) 

2C2 105 to 160 cm; brown (7.5TR 4/2) extremely gravelly sandy 


loam, yellowish brown (10YR 5/4) moist; massive; soft; few 


medium and coarse roots; 


65 percent rock fragments; neutral. 


Typic Ustifluvent, coarse-loamy, mixed, mesic (high water table). 


Pit location: 


Horizon 
Designation 





0 
AC 


Description 





T11S,R20W,SW1/4 NE1/4,Sec.4.San Francisco River. 


many very fine tubular pores; 


6 to 0 cm; litter layer of undecomposed leaves, twigs, etc. 


0 to 40 cm; light brownish gray (10YR 6/2) fine sand, brown 
(10YR 5/3) moist; single grain; loose; common very fine and 








Cl 


C2 


C3 


C4 


fine roots; many fine pores; laminar planes; strongly 
effervescent; moderately alkaline; abrupt smooth boundary. 


40 to 50 cm; light brownish gray (10YR 6/2) fine sandy 
loam, brown (10YR 5/3) moist; massive; loose; common very 
fine and fine roots; moderate fine and medium tubular 
pores; neutral; abrupt smooth boundary. 


50 to 127 cm; light brownish gray (10YR 6/2) fine sandy 
loam, brown (10YR 5/3) moist; single grain; loose, common 
very fine and fine, many coarse roots; many very fine and 
fine pores; slightly effervescent; moderately alkaline; 
abrupt smooth boundary. 


127 to 147 cm; light brownish gray (10YR 6/2) very gravelly 
sand, brown (10YR 5/3) moist; massive; slightly hard; few 
very fine pores; strongly effervescent; moderately alka- 
line; 40 percent rock fragments; abrupt smooth boundary. 


147 to 160 cm; light brownish gray (10YR 6/2) loamy very 
fine sand, brown (10YR 5/3) moist; single grain; loose; 
few fine pores; slightly effervescent; moderately alkaline. 


Aquic Ustifluvent, loamy-sketetal, mixed, mesic. 


Pit location: 


Horizon 
Designation 





AC 


2C1 


2C2 


3C3 


T11S,R20W,NE1/4 SE1/4,Sec.4,San Francisco River. 


Description 





0 to 7 cm; pale brown (10YR 6/3) cobbly loam fine sand, 
brown (10YR 4/3) moist; week granular structure; loose; few 
very fine roots; common very fine and fine pores; 15 per- 
cent rock fragments; slightly effervescent; moderately 
alkaline; abrupt smooth boundary. 


7 to 33 cm; pale brown (10 YR 6/3) extremely gravelly 
coarse sand, brown (10 YR 4/3) moist; single grain; loose; 
many very fine and fine, common medium roots, 80 percent 
rock fragments; slightly effervescent; moderately alkaline; 
abrupt smooth boundary. 


33 to 56 cm; pale brown (10YR 6/3) sand, brown (10YR 4/3) 
moist; single grain; loose; common fine and medium roots; 
slightly effervescent; moderately alkaline; abrupt wavy 
boundary. 


56 to 76 cm; pink (7.5R 7/4) silty clay loam; massive; 
hard; common very fine and many fine and medium roots; 
common medium tubular pores; strongly effervescent; 
moderately alkaline; abrupt smooth boundary. 


10 








7C7 


Typic Fluvaquent, lo 


Pit location: 


Horizon 
Designation 


Al 


2A2 


3C1 


4C2 


la 5 


76 to 88 cm; pale brown (10YR 6/3) extremely gravelly sand, 
brown (10YR 4/3) moist; single grain; loose; few very fine 
and fine roots; slightly effervescent; moderately alkaline; 
abrupt smooth boundary. 


88 to 107 cm; pink (10YR 7/4) silty clay loam, brown 
(10YR 4/3) moist; massive; hard; common fine and medium 
roots; few fine tubular pores; strongly effervescent; 
moderately alkaline; abrupt smooth boundary. 


107 to 135 cm; pale brown (10YR 6/3) extremely gravelly 
coarse sand, brown (10YR 4/3) moist; single grain; loose; 
many very fine and fine roots; common very fine and fine 
pores; 70 percent rock fragments; slightly effervescent; 
moderately alkaline; abrupt smooth boundary. 


135 to 160 cm; pale brown (10YR 6/3) extremely gravelly 
coarse sand, brown (10YR 4/3) moist; single grain; 70 
percent rock fragments; slightly efferescent; moderately 
alkaline; abrupt smooth boundary. 


amy-skeletal, mixed, mesic. 


T10S,R20W,SE1/4,SW1/4,Sec.8.San Francisco River. 


Description 


0 to 4 cm; light brownish gray (10YR 6/2) loam, brown 
(10YR 5/3) moist; weak fine platy parting to weak fine 
granular structure; slightly hard; many very fine roots; 
slightly effervescent; moderately alkaline; abrupt smooth 
boundary. 


4 to 13 cm; light brownish gray (10YR 6/2) fine sandy loam, 
brown (10YR 5/3) moist; weak fine granular structure; soft; 
many very fine roots; many very fine and fine pores; 
slightly effervescent; moderately alkaline; abrupt smooth 
boundary. 


13 to 70 cm; light brownish gray (10YR 6/2) extremely 
cobbly sand; brown (10YR 5/3) moist; single grain; loose; 
many fine and medium roots; strongly effervescent; 
moderately alkaline; abrupt wavy boundary. 


70 to 80 cm; light brownish gray (l10YR 6/2) sandy loam, 
brown (10YR 5/3) moist; massive; common fine and medium 
roots; strongly effervescent; moderately alkaline; abrupt 
wavy boundary. 


1] 














5C3 80 to 100 cm; light brownish gray (10YR 6/2) extremely 
cobbly coarse sand, brown (10YR 5/3) moist; single grain; 
loose; strongly effervescent; 80 percent rock fragments; 
moderately alkaline. 


We attempted to define the plant communities that occur within soil 
classes; weighted averages (WA) and reciproca’ averaging (RA) ordinations were 
constructed using ORDIFLEX (Gaugh 1977). RA is related conceptually to WA, 
except rather than applying preassigned species or plot weights, the method 
extracts relative weights for species and plots from the importance value 
matrix itself, using eigenvalue analysis. The result is a simultaneous 
ordination of species and plots that should correspond to some identifiable 
environmental gradient, in this case, gradations of hydric soils. We used the 
species ordering provided by RA to help construct a summary frequency table 
(Table 4) in which species are ordered along a soil moisture gradient from 
upper terrace to swale, with species frequency by soil indicated. This table 
allows quick determination of species that have indicator value with respect 
to differentiating among soils. "Obligate riparian" species, those species 
normally restricted to riparian or riparian-like habitats, were identified 
using criteria of Dick-Peddie and Hubbard (1977). Recent disturbance and 
subsequent succession on the floodplain can be inferred (Table 4) by the high 
frequencies of immature tree species (underlined numbers) such as Goodding 
willow (Salix goodingii) and Fremont cottonwood (Populus fremontii) in the 
short shrub layer of the sandbar and swale soils. Netleaf hackberry (Celtis 
reticulata) is common on the upper terrace soil as a tree, and it also has a 
large number of immature individuals on this same soil, as can be seen with 
frequencies of 13 in the tall shrub layer, 17 in the short shrub layer, and 
25, a high value, in the ground cover. 











The importance value weighted averages are given in Table 5. Importance 
value weighted averages with analysis of variance are given in Table 6, and 
presence/absence weighted average analysis is presented in Table 7. Tables 5, 
6, and 7 show that there is a strong correlation of high species frequency of 
occurrence index numbers with the drier upland soils, and the correlation 
continues with the lower index numbers and the progressively wetter soils. 


There is little difference between the “importance value" and "presence/ 
absence" analyses (Tables 6 and 7). Consequently, biomass information is 
probably unnecessary for this type of assessment. In all cases, the ANOVA's 
proved to have highly significant F ratios. Initially, we ran separate 
analyses for each river system (Gila and San Franciso), but the vegetation was 
so similar that pooling data to increase sample size seemed advisable. 


ORDIFLEX produced eight relatively discrete plant groupings that we 
classified as associations. These associations were derived independently of 
soil classes and serve as an independent assessment of soil class wetness 
(Table 8). The units also could be considered community types. There was a 
marshy spikerush (Eleocharis montevidensis) association and a barnyard grass 
(Echinichloa crusgalli) association. Both of these types were found on the 
swale soil, as would be expected from the frequency of occurrence index numbers 
of the species involved. Two associations were also found on the sandbar 
soil. These were the sandbar willow (Salix exigua) and the Goodding willow 
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Table 4. Frequency of occurrence of species on various soils along the Gila 


and San Francisco Rivers. 


Values are the percent occurrence of plants per 


soil class. Key to species code is included in appendix. 














Species Upper Lower 

code terrace terrace Sandbar Swale 

Tree jayer 

QUEM 50 

JUMO 47 4 

QUAR 42 . 

JUMA 38R 18R 

JUOS 27 4 

JUDE 27 2 

PIED 25 

FRYE 13R 8R 

CERE 10R 10R 

PRYI 7R 

QUGA 7 

ACNE 2R 8R 

MOMI 2R 

JUSC 2 

POFR 62R 22R 
SAGO 50R 7R 
PLWR 18R 2R 
POAC 8R 

POAN 4R 

ALOB 4R 

GLTR 2R 

Tall shrub layer 

RHTR 25R 4R 

FONE 13R,. 4R 

CERE 13R 12R 

HALA 13R 

PTAN 10R 

JUMO 7 

PRVI 5R 

FRYE 5 2R 
ALWR 2R 

GAWR 2R 

JUDE 2 

JUMA 2 6 
MOMI 2R 

PRGL 2R 

(Continued) 
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Table 4. (Continued) 

















Species Upper Lower 
code terrace terrace Sandbar Swale 
SAGO 12R 40R 
BAGL 10R 30R 
SAEX 6R 45R 
JUMA 6 
ACNE 4R 
POFR 4k 60R 
PLWR 4R 
TARA 2R 2R 
AMFR 2R 
SAIR 2R 13R 
Short shrub layer 
BRCA 50R 4R 
CERE 17R 12R 
JUMO 15 
JUMA 13 g 
VIAR 13R 2R 
HALA 10R 2R 
RHTR 10R 
FONE 10R 
GAWR 7R 
FRVE oR 2R 
QUAR 5 
JUDE 5 
PTAN 2R 
RHRA 2 
ALWR 2 
ARCA 2 
PIED 2 
PLWR 2R 
SAGO 4R 30R 23R 
POFR 40R 38R 
JUOS 2 2 
TARA 2R 3R 
ARLU 2 
ACNE 2R 
SAIR 5R 
FONE 2R 
(Continued) 
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Species 


code 


Table 4. 


(Continued) 





Upper 


terrace 


Lower 


terrace 


Sandbar 


Swale 





Ground cover 





STHI 
CERE 
QUAR 
MASP 
BRCA 
SENE 
HALA 
SPORO 
ERIGE 
BRICK 
GAWR 
VETH 
COCA 
BOCU 
LEPID 
PRY] 
MIGU 
MAYU 
VIAR 
LESQU 
JUMO 
PTAN 
PIED 
QUGA 
ERFL 
VIAL 
CONVO 
SOEL 
RHTR 
EVAL 
ARCA 
POA 
ERCA 
ERIOG 
FRANS 
GAMI 
FONE 
JUDE 
JUOS 
PEBA 


SIN ONS SS SY 
7A ZA 


yo) 
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Table 4. (Continued) 








Species Upper Lower 
code terrace terrace Sandbar Swale 





Ground cover (cont. ) 
POFE 
QUEM 
DAME 
BRMA 
ARLU 
CHENO 
TAOF 
HEDEO 
MEAL 16 17 
AMAR 14 35 l 
HOJU 
GAPA 
HEAN 
SAKA 
PAFL 
JUMA 
MESP 
BAHY 
PRGL 
SISYM 
CRTE 
GAURA 
MEOF 
AGAL 
SPCN 
SAEX 
MUHLE 
BAGL 
ARSP 
CLSE 
ERPO 
PANIC 
CONVO 
XASA 
EQAR + EQHY 7R 
POPE 67R 
SAGO 10R 23R 
AGSE 4R 
POMO 62R 
ECCR 42R 
JUTE 35R 
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Table 4. (Concluded) 





Species Upper Lower 
code terrace terrace Sandbar Swale 





Ground cover (cont.) 





ELMO 33R 
CECA 25R 
POAN 23R 
RACY 22R 
PORA 18 
TRIFO 15 
SCOL 10R 
PADI BR 
RONA 7R 
JUSA 7R 
JUTO 7R 
SEDUM 7 
TRRE 7 
CYPA 5R 
RAAQ 5R 
SCAL 3R 
PHPR 3 
TYDO 3R 
VEAN 3R 
RUCR 3R 
ACNE 3R 
TARA 3R 
TRLA 3 
ERME 2 
EVPI 2 
FRANS 2 
PLMA 2 
POAA 2 





aR" after a number indicates that the species is a riparian species as 
identified by Dick-Peddie and Hubbard (1977). 


Pan underline below a number indicates that the plant is an immature 
individual of a tree or shrub but is making a contribution to the vegetation 
of that layer. 
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Table 5. Importance value weighted averages by life-form for each soil class. 
Data are pooled for both rivers. 





Soil Mean Std.Err. N 

















Tree 
Upper terrace 4.44 0.141 40 
Lower terrace 2.02 0.084 50 
Sandbar 1.84 0.960 11 
Tall shrub 
Upper terrace 3.70 0.177 30 
Lower terrace 2.09 0.154 30 
Sandbar 1.49 0.070 40 
Swale 1.43 0.175 8 
Short shrub 
Upper terrace 3.93 0.127 36 
Lower terrace 2.44 0.243 25 
Sandbar 1.42 0.073 36 
Swale 1.71 0.030 34 
Ground cover 
Upper terrace 4.23 0.139 39 
Lower terrace 3.90 0.158 46 
Sandbar 3.71 0.158 40 
Swale 2.31 0.074 60 





associations. Again, the index numbers of the dominant species were consistent 
with these classifications. There was a considerable amount of California 
brickellia (Brickellia californica) on one of the sandbar soil replications, 
but it was not singled out as an association because it dominated only one 
replication. This same species was commonly a dominant of the shrub layer in 
lower terrace soil communities. Fremont cottonwood and Arizona sycamore 
(Platanus wrightii) were segregated as mature communities on the lower terrace 
soil of both rivers. Lanceleaf cottonwood (Populus acuminata), Arizona alder 
(Ainus oblongifolia), and velvet ash (Fraxinus velutina) were common 
subdominants in the Fremont cottonwood and Arizona sycamore associations. The 
upper terrace soil supported two associations, one on each river. An Arizona 
white oak (Quercus arizonica) association dominated this soil on the San 
Francisco River, whereas the upper terrace soil of the Gila River was dominated 
by an Emory oak (Quercus emoryi) association. In some stands, Arizona walnut 
(Juglans major) was common and almost a codominant species with the Arizona 
white oak. 
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Table 6. Importance value weighted averages with analysis of variance 

and Duncan's Multiple Range Test of weighted averages by soils and life-forms. 
Figures in parentheses indicate the number of plots involved. Data points 
joined by lines are not significantly different. 

















Tall Short Ground All life- 
Soil Tree shrub shrub cover forms 
Upper terrace 4.44(40) 3.70( 30) 3.93( 36) em 4.11(40) 
Lower terrace prea 2.09( 30) 2.44( 25) ages 2.70(50) 
Sandbar 1.84(11) poet genes 3.72(40) a 
Swale 1.42(8) 1.71(34) 2.32(60) 2.12(60) 
ANOVA 
F Value 83.83 55.92 23.12 44.25 114.55 
D.F. 2.00 3.00 3.00 3.00 3.00 
Significance .0001 .0001 .0001 .0001 .0001 
DISCUSSION 


We found a high correlation between hydric soils and wetland plants 
across moisture gradients on riparian ecosystems of the Gila and San Francisco 
Rivers of New Mexico. The three hydric soils all had importance value weighted 
average below 3 for all life forms combined (Table 6). Similar values were 
obtained for presence/absence weighted average analyses (Table 7). Vegetation 
on the nonhydric, upper terrace soils had average values over 4 by both methods 
of analysis. There is some question whether the lower terrace soils, with an 
importance value weighted average of 2.70, should be designated as hydric. 
Additional research is needed to determine the soil moisture relation of the 
lower terrace soils during the growing season, as they appear to be saturated 
for only brief periods, However, from our study, it would be valid to consider 
plant associations with lower index values as indicators of higher soil 
moisture, and associations with higher index values as indicators of drier 
soils. 
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Table 7. Presence/absence weighted averages with analysis of variance 

and Duncan's Multiple Range Test of weighted averages by soils and life-forms. 
Figures in parentheses indicate the number of plots involved. Data points 
joined by lines are not significantly different. 














Tal] Short Ground Ail life- 
Soil Tree shrub shrub cover forms 

Upper terrace 4.47(40) 3.77(30) 3.84( 36) 4. 33(39) 4.41(40) 
Lower terrace (1.89(50) 2.08( 30) 2.44( 25) ero 2.59( 40) 
Sandbar 1.82(11) any poem 3.67( 40) peices 
Swale 1.44(8) 1.76( 34) 2. 45( 60% 2.25(60) 
ANOVA 

F Value 22.08 60.49 71.99 53.16 123.81 

D.F. 2.00 3.00 3.00 3.00 3.00 

Significance .0001 .0001 .0001 .0001 .0001 





In addition to confirming a close correlation for the frequency of 
occurrence of species with low index numbers with hydric soils, the results of 
this work generally support the assignment of index numbers to New Mexico 
species by Reed (1986a). However, we encountered a number of species not on 
the New Mexico list that we believe should have an index number of 4 or lower. 
We assigned index numbers to those species as indicated in the Appendix. A 
number of those species are classified as “obligate riparian species" by 
Dick-Peddie and Hubbard (1987), and should be added to the State list of 
wetland plants for New Mexico. 


Dick-Peddie and Hubbard (1977) used the term “obligate riparian" for 
species that are limited to riparian or pseudo riparian habitats (ditch banks, 
foot of talus slopes, arroyos, and intermittent streams). In the Southwest, 
riparian habitats may not always have hydric soil conditions, as hydric soils 
are defined for this study. This is particularly true of soils associated 
with arroyos and intermittent canyon streams. Yet. the vegetation found 
associated with arroyos and intermittent streams is often composed of obligate 
riparian species, including some species, e.g., Fraxinus velutina and Tamarix 
ramosissima, that have been assigned hydric index numbers of 3 by Reed (1986a). 
The correct assignment of index number to all obligate riparian species wil] 
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Table 8. Riparian plant associations found on the four major soils of the 
Gila and San Francisco River floodplains. 








Soi] Plant community 
Upper terrace Arizona white oak association 
(Fluventic Ustochrept) ‘San Franciso River) 


Emory oak association 
(Gila River) 


lower terrace Fremont cottonwood association 
(Typic Ustifluvent) (both rivers) 


Arizona sycamore association 
(both rivers) 


Sandbar Sandbar willow association 
(Aquic Ustifluvent) (both rivers) 


Goodding willow association 
(both rivers) 


Swale Spikerush association 
(Typic Fluvenquent) (both rivers) 


Barnyard grass association 
(both rivers) 





facilitate identification and management of all western riparian ecosystems, 
wetland and otherwise. Our assigned numbers for these species appeared to be 
valid in the analysis. These species are often abundant enough to dominate or 
codominate their layers of riparian vegetation. We suggest that the following 
species, and index numbers, be added to the New Mexico list: tree--Acer 
negundo (2), Celtis reticulata (3), and Morus’ microphylla (3); 
shrub--Brickella spp. (4); forb--Erpilobium spp. (3), Parietaria floridana (2), 
Polygonum ramosissimum (3), and Stachys coccinea (2). 


























Results of the Gila and San Francisco Rivers study also suggest that some 
index numbers might warrant change. Fraxinus velutina (3) is often found on 
soils that are as wet as those found under Juglans major (2). On the other 
hand, in New Mexico, an index number of 2 is too low for Juglans major. We 
Suggest that the Juglans major index number be changed to 3. The same 
situation exists with Populus fremontii (2) and Salix qooddingii (1). These 
two obligate riparian species often occur as codominants and, in New Mexico, 
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Goodding willow is rarely found on the highly hydric soils typified by species 
carrying an index number of 1. We suggest that these two species carry the 
Same index number of 2. 


Forestiera neomexicana is an understory dominant in mature cottonwood 
riparian "gallery" forests and should carry the same index number (2) as 
cottonwoods. Helianthus annuus and Marrubium vulgare are ubiquitous species 
that colonize disturbed sites, and these two species have broad soil moisture 
tolerances. As a consequence, they have little indicator value for moist soil 


conditions, and should carry an index number no lower than 4. 
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APPENDIX 


Plants found on the Gila and San Franciso Rivers with frequency of occurrence 
(Reed 1986a) index numbers. 














Species Code Index No. 
Tree 
Acer negundo ACNE 2° 
Alnus oblongifolia ALOB 2 
Celtis reticulata CERE 
Fraxinus velutina FRYE 3 
Gleditsia triacanthos GLTR 3 
Juglans major JUMA 2 
Juniperus deppeena JUDE 53 
Juniperus monosperma JUMO 5 
Juniperus osteosperma JUOS 5 
Juniperus scopularum JUSC 5° 
Morus microphylla MOMI 3, 
Pinus edulis PIED 5 
Platanus wrightii PLWR 2 
Populus acuminata POAC 2 
Populus angustifolia POAN 2 
Populus fremontii POFR 2 
Prunus virens PRVI 38 
Quercus arizonica QUAR 5° 
Quercus emoryi QUEM 5 
Quercus gambelii QUGA 5° 
Salix gooddingii SAGO ] 
Shrub-tree 
Amorpha fruticosa AMFR 2 
Prosopis glandulosa PRGL 4° 
Ptelea angustifolia PTAN 4 
Tamarix ramosissima TARA 3 


(Continued) 
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Appendix. 


(Continued) 








Species Code Index No. 
shrub and vine 
Aloysia wrightii ALWR 4° 
Aplopappus laricifolius APLA 5° 
Artemisia carruthii ARCA 5° 
Artemisia ludoviciana ARLU 53 
Artemisia species ARTEM 5° 
Baccharis glutinosa BAGL . 
Brickellia californica BACA 4 
Brickellia fendleri BRFE 4° 
Brickellia species BRICK 4? 
Forestiera neomexicana FONE 4 
Gerrya wrightii GAWR 4? 
Rhus trilobata RHTR 4° 
Saiix exigue SAEX 1 
Salix irrorata SAIR 2 
Vitis arizonica VIAR 3 
Forb 
Ambrosia artemisifolia AMAR 4 
Bassia hyssopifolia BAHY 2 
Centaurium calycosum CECA 2 
Chenopodium species CHENO 5° 
Cicuta douglasii CIDO ] 
Cleome serrulata CLSE 3 
Conium maculatum COMA 2 
Convolvulus species CONVO 5° 
Conyza canadensis COCA 4. 
Croton texensis CRTE 25 
Cucurbita foetidissima CUFO 23 
Datura meteloides DAME 5 
Epilobium adenocaulon EPAD 3° 
Epilobium species EPILO 3° 
Erigeron divergens ERDI 5° 
Erigeron flagellaris ERFL 5° 
Erigeron species ERIGE 5° 
Eriogonum polycledon ERPO 5° 
Eriogonum species ERIOG 5. 
Erysimum capitatum ERCA a5 
Euphorbia albomarginata EVAL 4 
Evolvulus pilosus EVPI 53 
Franseria species FRANS 53 
Fragaria species FRAGA 3° 


(Continued) 








Appendix. (Continued) 








Species Code Index No. 
Galium microphy]]um GAMI 53 
Gaura parviflora GAPA 53 
Guara species GAURA 5. 
Hedeoma species HEDEO 5 
Helianthus annuus HEAN 3 
Helianthus species HELIA 53 
Lepidium medium LEME 53 
Lesquerella species LESQU x 
Macharanthera species MACHA 5 
Marrubium vulgare MAVU 3 
Melilotus alba MEAL 4 
Melilotus officinalis MEOF 4 
Mimulus guttatus MIGU 1, 
Mirabilis longiflora MILO m 
Parietaria floridana PAFL 2 
Penstemon barbatus PEBA 53 
Plantago major PLMA 2 
Polygonum persicaria POPE 2 
Polygonum ramosissimum PORA a 
Ranunculus aquatilis RAAQ hs 
Ranunculus cymbalaria RACY ] 
Rorippa nasturtium-aquaticum RONA ] 
Rumex crispus RUCR 2 
Salsola kali SAKA 4 
Senecio neomexicanus SENE 53 
Sedum species SEDUM 7 
Sisymbrium species SISYM 3 
Solanum elaeagnifolium SOEL ° 
Solanum rostratum SORO * 
Sphearalcea coccinea SPCO S 
Stachys coccinea STCO 2 
Taraxacum officinale TAOF 4. 
Trifolium lacerum TRLA 4 
Trifolium repens TRRE 4. 
Trifolium species TRIFO .. 
Urtica gracilenta URGR 3 
Veronica anagallis-aquatica VEAN 1. 
Verbascum thapsus VETH o 5 
Xanthium saccharatum XASA 4 


Grass and grasslike 
Agrostis alba AGAL 
Agrostis semiverticillata AGSE 





PRO PR 


(Continued) 
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Appendix. (Concluded) 








Species Code Index No. 
Bouteloua curtipendula BOCU 53 
Bromus japonicus BRJA 4 
Bromus marginatus BRMA 4 
Bromus rubens BRRU 43 
Cyperus parishii CYPA 2 
Echinochloa crusgal1i ECCR 2 
Equisetum arvense EQAR 2 
Equisetum hyemale EQHY 2 
Eragrostis mexicana ERME 4 
Hordeum jubatum HOJU 2 
Juncus saximontanus JUSA 2 
Juncus tenuis JUTE 2 
Juncus torreyi JUTO 2 
Muhlenbergia species MUHLE 23 
Paspalum distichum PADI l 
Panicum species PANIC 
Phleum pratense PHPR 4 
Poa annua POAN 3 
Poa species POA 53 
Polypogon monspeliensis POMO * 
Scripus acutus SCAC 1, 
Scirpus olneyi SCOL l 
Sitanion hystrix SIHY 5 
Sporobolu contractus SPCO 4 
Sporobolus cryptandrus SPCR .. 
Sporobolus species SPORO 4 
Typha domingensis TYDO ] 








“Plants not on Reed's (1986a) list. Frequency of occurrence index numbers 
assigned by the authors. 


29 


50272 - 10) 


























REPORT DOCUMENTATION }. REPORT NO. | 2. 3. Recipient's Accession No. 
PAGE | Biological Report 87(9) | 

4. Title and Subtitie 5. Report Date 
Soil-Vegetation Correlations on the Riparian Zones of the September 1987 
Gila and San Francisco Rivers in New Mexico 6. 

7. Author(s) — : . 8 Performing Organization Rept. No. 
W.A. Dick-Peddie, J.V. Hardesty, E. Muldavin & B. Sallach 

9. Performing Organization Name and Address _10. Project/Task/Work Unit No. 
Department of Biology _ 
New Mexico State University 11. Contract(C) or Grant(G) No. 
Las Cruces, NM 88003 (©) 14-16-0009-85-001 

(G) 
12. Sponsoring Organization Name and Address 13. Type of Report & Period Covered 


U.S. Department of the Interior 
Fish and Wildlife Service 





Research and Development 14. 
Washington, DC 20240 








15. Supplementary Notes 





16. Abstract (Limit: 200 words) 


This report documents the relationship between soils and vegetation in the riparian 
zone of the Gila and San Francisco Rivers in Southwestern New Mexico. Relationships 
were determined using the weighted average procedure developed by T.R. Wentworth and 
G.P. Johnson, North Carolina State University, for the U.S. Fish and Wildlife Service 
in 1986. Four soils, in order of proximity to the stream, were examined in this 
Study; swale, sandbar, lower terrace, and upper terrace. There was a positive 
correlation between hydric soils and wetland plants. Herbaceous ground cover was a 
less sensitive indicator of hydric soils than woody shrubs and trees. Soil moisture 
regimes were not well-defined, thus results should be considered tentative. 





17. Document Analysis a. Descriptors 


Riparian ecosystems 
Riparian vegetation 
Riparian soils 
Wetland ecology 


b. identifiers/Open-Ended Terms 


BEST COPY AVAILABLE 


New Mexico riparian ecosystems 
Southwestern U.S. riparian ecosystems 
New Mexico wetlands 


c. COSATI Fieid/Group 























18 Availability Statement 19. Security Class (This Report) 21. No. of Pages 
a Unclassified 29 
Release unlimited 20. Security Class (This Page) 22. Price 
Unclassified 
(See ANSI-—Z39.18) See Instructions on Reverse OPTIONAL FORM 272 (4-77) 


(Formerly NTIS—35) 
Department of Commerce 


% U.S. GOVERNMENT PRINTING OFFICE: 1987—777-216/65,329 REGION NO. 8 


Oo 











7% ae 
Hawauan isiands D 








i Headquarters, Research and Development, 


Washington, DC 


@ Locations of Regional Offices 


REGION 1 

Regional Director 

U.S. Fish and Wildlife Service 

Llovd Five Hundred Building, Suite 1692 
500 N.E. Multnomah Street 

Portiand, Oregon 97232 


REGION 4 

Regional Director 

U.S. Fish and Wildlife Service 
Richard B. Russell Building 
75 Spring Street, S.W. 
Atlanta, Georgia 30303 














REGION 2 

Regional Director 

U.S. Fish and Wildlife Service 
P.O. Box 1306 

Albuquerque, New Mexico 87103 


REGION 5 

Regional Director 

U.S. Fish and Wildlife Service 

One Gateway Center 

Newton Corner, Massachusetts 02158 


REGION 7 

Regional Director 

U.S. Fish and Wildlife Service 
1011 E. Tudor Road 
Anchorage, Alaska 99503 


BEST COPY AVAILABLE 











Puerto Rico and 
.-@* 
C7. 
Virgin islands 











REGION 3 


Regional Director 

U.S. Fish and Wildlife Service 
Federal Building, Fort Snelling 
Twin Cities, Minnesota 55111 


REGION 6 

Regional Director 

U.S. Fish and Wildlife Service 
P.O. Box 25486 

Denver Federal Center 
Denver, Colorado 80225 


3\ 











Preserve Our Natural Resources 


U.S. 
FISH @ WILDLIFE 
/ERVICE 


DEPARTMENT OF THE INTERIOR |[/, 
U.S. FISH AND WILDLIFE SERVICE 


As the Nation's principal conservation agency, the Department of the Interior has respon- 
sibility for most of cur nationally owned public lands and natural resources. This includes 
fostering the wisest use of our land and water resources, protecting our fish and wildlife, 
preserving the environmental and cultural values of our national! parks and historical places, 
and providing for the enjoyment of life through outdoor recreation. The Department as- 
sesses Our energy and mineral resources and works to assure that their development is in 
the best interests of ali our people. The Department also has a major responsibility for 
American indian reservation communities and for people who live in island territories under 
U.S. administration. 












BEST COPY AVAILABLE 


%¢ 











FILMED 
04-03-89 








